Analysis of a time discretization scheme for 
a nonstandard viscous Cahn-Hilliard system^ 

Pierluigi Colli'^^) 

^^ ' e-mail: pierluigi.colli@unipv.it 

(N ! Gianni Gilardi^^^ 

^ I e-mail: gianni . gilardiSunipv . it 



< 



> 



% 



Pavel Krejci^^^ 

e-mail: krejci@math.cas.cz 

Paolo Podio-Guidugli'^'^^ 

e-mail: ppg@unironia2 . it 



< 

^ ; Jiirgen Sprekels^^^ 

d ' e-mail: sprekels@wias-berlin.de 



^ ' Dipartimento di Matematica "F. Casorati", Universita di Pavia 
Via Ferrata 1, 27100 Pavia, Italy 



\— ^ I "^ Institute of Mathematics, Academy of Sciences of the Czech Republic 

l/^ ' Zitna 25, CZ - 115 67 Praha 1, Czech Repubhc 

• . ^ ' Dipartimento di Ingegneria Civile e Ingegneria Informatica 

f^ , Universita di Roma TorVergata, Via Politecnico 1, 00133 Roma, Italy 

m 



' ' Weierstrafi-Institut fiir Angewandte Analysis und Stochastik 
MohrenstraBe 39, 10117 Berlin, Germany 

Abstract 

In this paper we propose a time discretization of a system of two parabolic equa- 
tions describing diffusion-driven atom rearrangement in crystalline matter. The 
equations express the balances of microforces and microenergy; the two phase fields 
are the order parameter and the chemical potential. The initial and boundary- value 
problem for the evolutionary system is known to be well posed. Convergence of the 
discrete scheme to the solution of the continuous problem is proved by a careful 
development of uniform estimates, by weak compactness and a suitable treatment 
of nonlinearities. Moreover, for the difference of discrete and continuous solutions 
we prove an error estimate of order one with respect to the time step. 
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2 Time discretization of a nonstandard viscous Cahn-Hilliard system 

1 Introduction 

This paper deals with a time discretization of a PDE system arising from a mechanical 
model for phase segregation by atom rearrangement on a lattice. The model was proposed 
by one of us in [29]; the resulting system has been analyzed in a recent and intensive 
research work by four of the present authors (see, in particular, [16] and [17], both for 
well-posedness results and for a detailed presentation of the model). 

The initial and boundary value problem we consider consists in looking for two fields, 
the chemical potential fi > and the order parameter p G (0, 1), solving 



l + 2g{p))dtfi + fidtg{p)-Afi = 


infix (0,T), 


(1.1) 


dtp-Ap + f'{p) = fig'{p) 


infix (0,T), 


(1.2) 


dup. = dup = 


onT X (0,r), 


(1.3) 


/w|i=o = Ato and p\t=o = Po 


in fi; 


(1.4) 



here fi denotes a bounded domain of M^ with conveniently smooth boundary F, T > 
stands for some final time, and di, denotes differentiation in the direction of the outward 
normal u. 

Problem (1.1)-(1.4) is parameterized by two nonlinear scalar-valued functions, g and 
/, which enter into the definition of the system's free energy: 

i, = z^(p, Vp, /") = - (^ + 9{p}j^^ + fip) + ^IVpp. (1.5) 

We point out that in (1.1)-(1.5) all physical constants have been set equal to 1. We 
also note that the last two terms in (1.5) favor phase segregation, the former because 
it introduces local energy minima, the latter because it penalizes spatial changes of the 
order parameter. For g, one can take any smooth function, provided it is nonnegative in 
the physically admissible domain: 

g{p) >0 for all pE (0,1); (1.6) 

accordingly, the coefficient 1/2 of /i in (1.5) should be regarded as a prescribed material 
bound. As to the possibly multi-well potential /, we take it to be the sum of two functions: 

/(p) = /i(p) + /2(p); 

the one, /i, is convex over (0, 1), and such that its derivative /( (and possibly also /i) is 
singular at the endpoints and 1 (cf. (2.3)); the other is required to be smooth over the 
entire interval [0, 1], but not to have any convexity property, so that in equation (1.2) /g 
may serve as a non-monotone perturbation of the increasing function /(. 

As to the parameter functions, in [16] the choice made for g was: 

9ip)=P, (1-7) 

while the assumptions on / were compatible with choosing a double-well potential: 

/(p) = ai {p ln(p) + (1 - p) ln(l - p)} + ^2 p (1 - p) + asP, (1-8) 



A;u„+i = 


in Q, 


(1.9) 


/in^'(Pn) 


in fi, 


(1.10) 





on r, 


(1.11) 
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for some non-negative constants Oi, a2, as- Note that, if 03 is taken null, then, according 
to whether or not 2ai > 02, / turns out to be convex in the whole of [0, 1] or it exhibits 
two wells, with a local maximum at p = 1/2; moreover, for a^ > 0, the combined function: 

-9ip)l^ + fip) (a part of ip) 

shows one global minimum in all cases, and it depends on the sign of [a^ — /i) which 
minimum actually occurs. On the other hand, the framework of paper [17] allows for 
much more general choices of g and /, as well as for nonlinear diffusion of /i. Existence 
and uniqueness results were proved in both [16] and [17], with different approaches. Here, 
we take inspiration from arguments developed either in the one or in the other of those 
papers. 

We introduce a time discretization of system (1.1)-(1.4) which is implicit with respect 
to the principal terms and tries to handle very carefully the nonlinearities. Namely, we 
address the recursive sequence of the elliptic problems: 

(1 + 27„) ShfJ^n + Pn+l Shin 

6hPn - Ap„+i + f'{pn+l) 

duPn+l = duPn+l 

ioT n = 0,1, . . . ,N, where h = T/N is the time step, jn '■= g(Pn) and, for any (A^ + l)-ple 

Zq, z\, . . . , Zn, we let 

5hZn ■= [Zn+l - Zn)/h for 72 = 0, . . . , iV - 1. 

After showing the existence of a discrete solution at any step, we carry out a number 
of uniform estimates on the time-discrete solution which allow us to prove convergence 
to the unique solution (p, p) of the continuous problem (1.1)-(1.4), as h tends to (or, 
equivalently, N goes to +00). Then, we estimate certain norms of the difference between 
the piecewise-linear-in-time interpolants of the discrete solutions and the continuous so- 
lution: more precisely, the first error estimate we prove is of order /i^/^; the second, which 
holds under stronger regularity assumptions on the initial data, is of order h. 

We regard our results as a cornerstone in the construction of a time-and-space dis- 
cretization of problem (1.1)-(1.4). With reference to such a complete discretization of 
Cahn-Hilliard and viscous Cahn-Hilliard systems, we quote papers [1, 2, 3, 4, 5, 6, 7, 8, 
22, 21, 23]. Some recent efforts can be found in the literature with the aim of analyzing 
other classes of phase transition problems, either to show existence via time discretiza- 
tion [9, 14, 15, 19, 20, 27, 30, 35, 36] or to prove numerical results such as special conver- 
gence properties, stability or error estimates [11, 12, 13, 18, 25, 28, 31, 33, 34] (cf. also 
[26] for a recent review on phase-field models). We dare say that our contribution goes 
deeply into the structure of the mathematical problem, because, as is not the case for 
many other similar investigations, we succeed in showing a linear order of convergence. 

Our paper is organized as follows. In the next section, we list and discuss our as- 
sumptions, formulate the continuous and discrete problems precisely, and state our main 
results. Section 3 is devoted to proving that there is a discrete solution. The convergence 
result is proved in the long and articulate Section 4. Finally, the last two Sections 5 and 6 
contain detailed proofs of the two error estimates. 
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2 Main results 

In this section, we describe the mathematical problem under investigation, introduce the 
time discretization scheme, make our assumptions precise, and state our results. 

First of all, we assume ^2 to be a bounded connected open set in M^ with smooth 
boundary F. For convenience, we set 

V := H^{Vi), H := ^^(fi), and W := {v G H^{n) : d^v = 0, on F}, (2.1) 

and we endow these spaces with their standard norms, for which we use a self-explanatory 
notation like || • \\v- The notation || ■ \\p {1 < p < +oo) stands for the standard L^-norm 
in L^{Q); for short, we sometimes do not distinguish between a space (or its norm) and a 
power thereof. 

As to the parameter functions / and g, we assume that 

/ = /i + /2, where (2.2) 

/i : (0, 1) — 7- [0, +oo) is a convex C^ function satisfying 

lim /{(r) = — oo , and lim/((r) = +cxd, (2.3) 

/s : [0, 1] ^ M is of class C^; (2.4) 

(7 : [0, 1] — 7- M is of class C^ and nonnegative. (2.5) 

For the initial data, we require that 

/io e 1/ n L°°{Q) and /io > a.e. in Q; (2.6) 

PoEWcC^m and inf po > 0, suppo < 1- (2.7) 

We stress that conditions (2.7) actually imply that po is 1/2-Holder continuous: indeed, 
as n is a three-dimensional domain, W is continuosly embedded in C^'^^'^{Q.). As a 
consequence, also /(po) and /'(po) are 1/2-Holder continuous, since / and /' are smooth 
in (0, 1). On the other hand, we point out that in the sequel we will mostly exploit the 
compactness of the embedding W C C°(r2); Holder continuity will play no role. 

As recalled in the Introduction, in papers [16] and [17] two versions of problem (1.1)- 
(1.4) were solved over an arbitrary time interval [0,T] in a rather strong sense, because 
the solution pairs (/i, p) were required to satisfy 

fieH\0,T;H)r]L\0,T;W), (2.8) 

p G W^'°°{0, T; H) n iJ^(0, T; V) n L°°(0, T; W), (2.9) 

p > a.e. in Q, (2.10) 

0<p<l a.e. in Q and f'{p) E L'^{0,T; H). (2.11) 

Note that the boundary conditions (1.3) follow from (2.8)-(2.9), due to the definition 
of W in (2.1). Accordingly, the solutions to the problems of type (1.1)-(1.4) studied in 
[16] and [17] were pairs (pt, p) satisfying, in addition to (2.8)-(2.11), the system 

{l + 2g{p))dtp + pdtgip)-Ap = a.e. in Q, (2.12) 

dtp - Ap + f{p) = fig'ip) a.e. in Q, (2.13) 

p(0) = Po and p(0) = po a.e. in fi. (2-14) 

Some of the results proved in the quoted papers are summarized in the following theorem. 
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Theorem 2.1. Let assumptions (2.2)-(2.7) hold. Then, there exists a unique pair {fi,p) 
satisfying (2.8)-(2.11) and solving problem (2.12)-(2.14). Moreover, fi G L°°{Q), and 
there exist p^, p* G (0, 1) such that P:^ < p < p* a.e. in Q. 

The main aim of the present paper is to show that, given the time-discretization 
scheme introduced here below, the discrete solution converges to the solution (/i, p) as the 
time step h tends to zero. 



Notation 2.2. Assume that A^ is a positive integer, and let Z be any normed space. We 
define 6h : Z^+^ -^ Z^ as follows: 

for z = {zq, zi, . . . , zn) e Z^^^ and w = (wq, . . . , wn-i) E Z^, 

5hZn = w means that Wn=— — ; for n = 0, . . . , A^ — 1. (2-15) 

h 

We can also iterate the discretization procedure, and define, e.g., 

Oh^n-= 7 = p forn = 0,...,A^-2. (2.16) 

Next, by setting h := T/N (without stressing the dependence of h on A^) and /„ := 
((n — l)h,nh) for n = 1,. . . ,N, we introduce the interpolation maps from Z^"*"^ into 
either L°°(0,T; Z) or W^'°^{0,T; Z) as follows: for z = {zq^Zi, . . . , Zn) e Z^+\ we set 

Zh,Zf^eL°°{0,T;Z) and 2;, G H^^'~(0, T; Z) , (2.17) 

Zh{t) = Zn and z_}SS) = -^n-i for a. a. t E In, n = 1, . . . ,N, (2.18) 

Zh{0) = Zq and dtZh{t) = ShZn-i for a.a. t G /„, n = 1, . . . , A^. (2.19) 

These maps yield the backward/forward piecewise-constant and piecewise-linear inter- 
polants of the discrete vectors. We obviously have: 

\\Zh\\L^(0,T;Z) = max \\Zn\\z, |khlU°°(0,T;Z) = HiaX \\Zn\\z, (2.20) 

n=l,...,N n=0,...,N—l 

N N~l 

11-^ ||2 _/,V^|^ |2 11^ ||2 _/,\^|~ |2 (O 0T\ 

\\'^h\\L'2(0,T;Z) — ''' /_^\'^n\z ^ ll^/i IIl2(0,T;Z) ~ '' _^ F" Iz • {Z.Zl) 

n=l n=0 

Moreover, as ^h{t) is a convex combination of Zn-i and Zn for t G /„, we also have 

\\zh\\L^{o,T;Z) = ^max^max{||2„_i||z, \\zn\\z} = max{||2;o|U, \\zh\\L°°{o,T;Z)}, (2.22) 

N 

lk/i|lL2(0,T;Z) ^ ^2^(lkn-l|lz + lkn|U) < ^lko|U + 2|| 2;h||L2(0,T;Z) • (2.23) 

n=l 

Finally, by a direct computation, it is straightforward to prove that 

\\Zh — Zh\\L°°(0,T;Z) = max \\Zn+l — Zn\\z = hWdtZhWLoo/Q^rp.^z), (2.24) 

n.=0,...,A' — 1 

h^~^ h^ ^ 

\\Zh — Zh\\L'i{Q^T;Z) — TT _^ll^n+l — Zn\\z = -^ \\^tZh\\L'^{Q^T;Z)-i (2.25) 

n=0 
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and that the same identities hold for the difference z_j^ — ^h- 

At this point, we can write the discrete scheme presented in the Introduction in a 
precise form. For any positive integer N, we look for two vectors {fin)n=o ^^^ {pn)n=o 
satisfying the following conditions: 

i) the first components /io and po coincide with the initial data; 

a) for n = 0, . . . ,N — 1, we have that 

/i„+i , pn+i e W, Hn+i > and < pn+i < 1 in fi, /'(p„+i) G H; (2.26) 

in) if {'yn)n=o ^^ ^^^ vector whose components are 7„ := g{pn), there hold 

(1 + 27„) (5/,/i„ + yu„+i Shin - Ayu„+i = 0, (2.27) 

5hpn - Ap„,+ i + /'(p„,+ l) = Png'ipn), (2.28) 

forn = 0,...,N -1. 

Also in this case, the homogenous Neumann boundary conditions are implicit in the 
regularity requirements (see (2.26) and (2.1)). 

Clearly, the "true" problem consists in finding {pn+i, Pn+i) once (p„, pn) is given. Here 
is our result in this direction. 

Theorem 2.3. Assume (2.2) -(2.7). Then, there exists h^ > such that, for h < Hq 
and n = 0, . . . , A^ — 1, problem (2.27)-(2.28) has a unique solution {pn+i, Pn+i) satisfing 
(2.26). 

Our next results concern firstly convergence of interpolants for vectors (p„,) and (/i„) 
to the solution {p,p) to problem (2.12)-(2.14), then error estimates. We point out that, 
for simplicity, the convergence theorem here below is not stated in a precise form: the 
topological setting will be specified later, by means of relations (4.34)-(4.38). 

Theorem 2.4. Assume (2.2)-(2.7). Then, in accord with Notation 2.2, the sequences of 
interpolants for the discrete solutions given by Theorem 2. 3 converge to the solution (/i, p) 
given by Theorem 2.1 as h tends to 0, in a suitable topology. 

Theorem 2.5. In addition to (2.2)-(2.7), assume that 

po G H^in). (2.29) 

Then, for sufficiently small h > 0, the following error estimate holds: 

\\Ph — P\\m{0,T;H)nL°°{0,T;V) + WPh — P\\L°°(0,T;H)nL2{0,T;V) < ck ' , (2.30) 

where c depends only on the structural assumptions and the data. 
Theorem 2.6. In addition to (2.2)-(2.7), assume (2.29) and 

Po e W. (2.31) 

Then, for sufficiently small h > 0, the following error estimate holds: 

\\Ph — P\\H'^{0,T;H)nL'^{0,T;V) + \\Ph " P\\L°°{0,T;H)nL^{0,T;V) < C /i, (2.32) 

where c depends only on the structural assumptions and the data. 
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Remark 2.7. It is easy to see that our assumptions (2.2)-(2.7) ensure that both /'(po) 
and /iofi''(po) belong to V. It follows that (2.29) is equivalent to 

-Apo + /'(po)-W(po)eV^. (2.33) 

We also notice that the assumptions (2.29) and (2.31) ensure further regularity for the 
solution {fi, p) to the continuous problem (see the forthcoming Remark 6.1). 

We prove Theorem 2.3 in Section 3 and Theorem 2.4 in Section 4; the last two sections 
are devoted to proving, respectively, Theorems 2.5 and 2.6. 

Throughout the paper, we account for the well-known embeddings V C L'^{Q) (1 < 
g < 6) and W C C'^^Q), and for the related Sobolev inequalities: 

\\v\\q < C||f llv and ||t'||oo < C'IK'IIvk, (2.34) 

ioT V & V and v G W, respectively, where C depends on Q only, since sharpness is not 
needed. We remark that these embeddings are compact. In particular, the following 
compactness inequality holds: 

II^IU < "^llVf IIh + Co-||t'||_f/, for every v E V and a > 0; (2.35) 

in (2.35), C„ is a constant that depends only on Q and a. Furthermore, we make repeated 
use of Holder inequality, of the following elementary identity: 

[a-b)a= -a^ --b"^ + -{a-bf, for every a, 6 e M, (2.36) 

and of Young's inequality 

ab < aa^ H b"^, for every a, 6 > and a > 0. (2.37) 

4cr 

Moreover, we use the discrete Gronwall lemma in the following form (see, e.g., [24, 
Prop. 2.2.1]): if (oq, . . . , a^v) e [0, +00)^+^ and (61, . . . , 6iv) e [0, +cx))^ satisfy 

771—1 



0777 < «o + y^ bndn for 771 = 1, . . . , A^, then 

77 = 1 

■m— 1 

^777 < Oo exp ( 2. ^n) for m = 1, . . . , A^. (2.38) 



71=1 

Finally, throughout the paper we use a small-case italic c for a number of different 
constants that may only depend on fi, the final time T, the shape of /, the properties 
of the data involved in the statements at hand; those constants we need to refer to are 
always denoted by capital letters, just like C in (2.34). Moreover, a notation like c^ signals 
a constant that depends also on the parameter a. The reader should keep in mind that 
the meaning of c and Cg. might change from line to line and even in the same chain of 
inequalities and that their values never depend on the time step h. 
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3 Existence 

In this section, we prove Theorem 2.3. We argue inductively with respect to n, i.e., by 
assuming that a pair (/i„, p„) satisfying (2.26) with n in place of n+1 is given, we prove that 
problem (2.27)-(2.28) has a unique solution (yU„+i,p„+i) satisfying (2.26). More precisely, 
as is going to be clear from the proof, we need less regularity for //„, e.g., //„ G V . In 
particular, our assumptions on /io are sufficient to start. We rewrite (2.27)-(2.28) in the 
form 

(l + 7„ + 7„+i)/i„+i - hAfin+i = (l + 27„)/in, (3.1) 

Pn+l - hApn+1 + hf'{pn+l) = Pn + hpng'{Pn), (3.2) 

and solve (3.2) first for p„+i (so that jn+i is also known), then (3.1). In order to solve 
both problems, it is expedient to replace each equation by a minimum problem, at least 
for h small enough. We consider the functionals: 

Ji : y — i> M and J2 : V -^ (—00, +00], defined by, respectively, 

^1(^0 •= o / |Vtf + ^ / (1 + 7n + In+iy - / (1 + 27n)/in V and (3.3) 

Uv) := ^ / |Vtf + ]- I v^ + h f f{v) - f{pn + hfi^g'ip^)) v. (3.4) 

-^ Jn ^ Jn Jn Jn 

In (3.4), we have f = fi + /2, where /2 is any smooth extension of /2 to the whole 
of M and /i is the unique convex and lower semicontinuous extension of /i that satisfies 
/(r) = +00 if r ^ [0, 1]. By the way, it is understood that the corresponding integral that 
appears in (3.4) is infinite if f{v) does not belong to L^{Q). Therefore, both functionals 
are well-defined and proper whenever pn & V and pn & W (and this implies boundedness 
of g{pn) and g'{pn))- Moreover, in view of the above remarks, Ji is continuous, and J2 
lower semicontinuous, on V. 

Now, we observe that equations (3.1) and (3.2), when complemented by the regularity 
requirements in (2.26) (which yield the homogeneous Neumann boundary conditions), 
are the strong forms of the Euler-Lagrange variational equations for the stationary points 
of Ji and J2, respectively. More precisely, the strong form (3.1) follows from the varia- 
tional formulation thanks to the regularity theory of elliptic equations. As far as (3.2) 
is concerned, the function /' should be replaced - in principle, at least - by the sum 
dfi + /2, where dfi is the sub differential of /i. However, once an L^(r2)-estimate is ob- 
tained for the sub differential (and standard arguments of the theory of maximal monotone 
operators (see, e.g., [10]), easily yield such an estimate), the variational Euler-Lagrange 
equation can be written exactly in the form (3.2), because dfi is single- valued due to our 
assumptions on /i (see (2.3), in particular). Consequently, existence and uniqueness of 
the solution [pn, pn) follow if the functionals (3.3) and (3.4) are convex, so that each of 
the correponding minimum problems has a unique solution. This is granted for the first 
problem: indeed, Ji is strictly convex and coercive, because g is nonnegative. The same 
holds for J2, provided that the second derivative of function r i-> r^/2 -|- hf2{r) is strictly 
positive on [0, 1], which is the case if /isup l/2'l < 1. 

It remains to prove that Pn+i > 0. To this end, we multiply (3.1) by —p~^i, where 
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v~ = max{— u, 0} denotes the negative part of v, and integrate over Q. We obtain: 

/ (1 + gipn) + g{pn+l)) l/^n+lP + ^ / |V/i;;+iP = - / (1 + 2g{pn))fInPn+l < 0) 

Jq Jn Jn 

because both g and /i„ are nonnegative. This imphes that Pn+i = 0? ^^-iid hence that 
/i„+i > 0. 

4 Convergence 

In this section, we prove Theorem 2.4. For convenience, we introduce one more vector, 
{^n)n=oy ^^^ recall the definition of (7n)^=o- 

Cn := f[ipn) and 7„ := g{pn) forn = 0, . . . , iV. (4.1) 

Later on, we also use the interpolants of these vectors according to Notation 2.2. Our 
argument uses compactness and monotonicity methods. 

First a priori estimate. We multiply (3.1) by /in+i and integrate over ^2. By account- 
ing for the elementary identity (2.36), we obtain 

- / /i^+i " 2 / ^" ^ 2 / '^''^^ ~ ^'''^ ^ ^ / l^^"+il^ 

+ / {inPl+l + 7n+l/i^+l - 27„/i„/i„+i) = 0. 

As 7n/in+i + 7n+i/i^+i " 27„/i„/i„+i = 7„+i/i2^^ _ ^^^2 _^ ^^(^^_^^ - jj^f ^ we derive that 

+ /l / |V;U„+i|2 = 0. 

On summing over n = 0, . . . , m — 1 with 1 < ?72 < A^, we conclude that 

m~l „ , m— 1 



/* -| / ' (■ X /i ^ lit _1_ /i 

/ (- + 7m) /i^ + ^' 5^ / (o + ^") l'^'^^"!' + ^ Zl / 1^^ 

o+7o)Ato for m= 1,...,7V. 



l2 
n+1 



n=0 "^^^ n=0 

1 

As (7 is nonnegative and hence 7j > 0, this implies that ||/im||H < c for m = 1, . . . , iV. 
Thus, the above estimate also yields 

N~l N 

max Wl^mWn + h'^y^JShPnWH + hy^WiinWv <c. (4.2) 

n=0 n=l 

In terms of the interpolants introduced in Notation 2.2, with the help of /xq ^ V, (2.20)- 
(2.21), and (2.24)-(2.25) we have that 

11— ||2 I II ||2 

\\f^h\\L°°{0,T;H)nL^iO,T;V) ' \\t!Lh\\L°°{0,T;H)nL^{0,T;V) 

+ \\P'h\\L^[o^T;H)nL'2{0,T;V) + M'^tPhW L^(o,T;H) ^ C. (4.3) 
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Second a priori estimate. In (3.2), we move p„ to the left-hand side. Then, we 
multiply by pn+i — Pn and integrate over Q. We obtain 

Pn+l - Pn? + 77/ \^Pn+l? " 7T / \^ Pn? + 7T / I^Pn+l " Vp„|^ 



2 1 I / it,-r i- 1 f) I I / 't I Q 

Jn '^ Jn ^ Jn 



+ h I f'{pn+l)ipn+l - pn) = h / Pn9' ipn){pn+l - Pn) ■ (4.4) 



Now, we consider the last integral on the left-hand side of (4.4). We split f = f[ + /2 
and use the convexity assumption of /i and boundedness for /g. We get 

/'(Pn+l)(Pn+l - Pn)> / (/l(pn,+l) " /l(Pn,)) " C / |p„+i - p„ | . 

Since also g' is bounded, we infer from (4.4) that 



/ \pn+l - pn? + 2 / l^/^"+ll^ ~ 2 / 1^^"'^ ^ 2 / 1^^"+^ ^ ^P"" 
o \i t/iZ i/ii t/'Z 



|2 



+ /i / (/l(Pn+l) - /l(Pn)) < c/l / (1 + /i„)|p„+i - p„| 
< 2 / '^"+1 -pnP + c/i^ {I + pD <- \pn+l - Pn? + Ch^, 

the last inequality being due to (4.3). By dividing by h, summing over n = 0, . . . ,m — 1, 
and owing to the obvious inequality mh < c, we conclude that 

hT / l^hPnP + - / |Vp„,+i|2 + h^y^ \5hiyPn)? + / /l(p„^+l) < C (4.5) 

n=o-^^ ^>/n ^^o^n Jn 

for 7T7, = 0, . . . , iV — 1. As the term involving the difference quotient 6hPN-i is missing in 
the first sum since ?Ti < A^ — 1, we estimate it directly. We multiply (2.28), written for 
n = A^ — 1, by h6hPN-i and integrate over fl. We have 

h / \ShpN-l? + / (VpAT - VpN-l) ■ Vpjv + / f[{pN){pN - pN-l) = h (f)6hpN-l, 



where we have set (p := Pn-i g'{PN-i) — f2{PN-i)- Owing to the elementary identity (2.36) 
and to the convexity of /i as before, we have 

11 f 

h\\ShpN-i\\H + 2 W^PnWh + 2 W^P^ ~ '^Piv-illn + / /i(P7v) 

t/ i z 

< - ||Vpjv-i||h + / /i(pjv-i) + h\\c\>\\n W^Pn-iWh 
^ Jn 

< 2 I|Vpjv-iIIh + / h{pN-i) + -Wh + i:¥hPN^i\\H- 

Now, we observe that the first two terms of the last line are bounded by (4.5) written with 
m = N — 2 and that (p is estimated in H thanks to (4.2) and our assumptions of g and f2- 
Moreover, the last term of the first line can be ignored since /i is nonnegative. Hence, we 
get the desired bound for the first term. At this point, we can easily derive an estimate 
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for IIptoII// for m = 1, . . . , A^. By using the obvious identity Pm = Po + h X]™=o ^hPn and 
the euchdean Schwarz and Young inequahties, we see that 

m—l , m—1 

WPtuWh < WPoWh + h^\\6hPn\\H < C+ -\m+ ^WShPnWHj < c. 
n=0 n=0 

Hence, by recaUing (4.5) and our last estimates, we conclude that 

7V-1 Af-2 

hY.\\^hPn\\l+ max \\pm\\v + h^Y.\\^h{^Pn)\\l<c. (4.6) 

•'— ' m=l,...,N •'— ' 

n=0 n=0 

In terms of the interpolants, (4.6) reads (thanks also to po ^ ^ and to (2.20)-(2.21) and 
(2.24)) 

ll^tP/i|lL2(0,T;H) + \\Ph\\L°°{0,T;V) + IIP,jIIl°°(0,T;V) 

+ l|Pft.||L°°(0,r;y) + ^ II ^t^Pft II L2(o,T-h;/f) < C. (4.7) 

Third a priori estimate. We come back to (2.28) and rewrite it as (recall (4.1)) 

-Ap„,+i + ^„+i = -6hPn + Png'iPn) " /2(Pn+l)- 

Hence, a standard argument (multiplying by — Ap^+i and by ^„+i) shows that the follow- 
ing estimate holds true 

||Ap„+i||h + ||^„+i||h < c\\-6hPn + Png'iPn) " /2(Pn+l)||H • 

Thus, we infer that 

II Apn+illi + ll^n+ill^ < c{\\6h.Pn\\H + llPnll^ + l) for n = 0, . . . , A^ - 1. (4.8) 

Moreover, by using the regularity theory of elliptic equations, we deduce that 

||Pn+l||^ + Un+l\\l < c{\\Pn+l\\v + W^hpnWH + Wf^nlln + l) • (4-9) 

Now, we multiply (4.8) by h and sum over n = 0, . . . ,m — 1. By accounting for (4.2) 
and (4.6), we conclude that 

N-l N-1 

hJ2\\P-+^\\w + hJ^Un+iWl < c. (4.10) 

n=0 n=0 

In terms of the interpolants, (4.10) yields (by accounting for po G W) 

IIP/illL2(o,r;Vl/) + \\pJ\l^O,T;W) + IIP/i|Il2(0,T;H/) ^ ^ (4-11) 

besides an estimate for, e.g., ||^^|| in L'^(0,T; H). 

Fourth a priori estimate. We write (2.28) with n + 1 in place of n and take the 
difference between the obtained equality and (2.28) itself. Then we multiply this difference 
by ShPn+i and integrate over fi. We have 



/ (ShPn+l - ShPn)Shpn+l + / (Vpn+2 - Vpn+l) " 5hV pn+1 

Jn Jn 



+ / {f'{Pn+2) - f'{pn+l))ShPn+l = / {P'n+ig' (pn+l) - PuQ' {pn))ShPn+l- (4.12) 

Jn Jn 
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By accounting for the elementary identity (2.36), we get 

/ (ShPn+l - ShPn)ShPn+l = TT / l^/iPn+l^ ~ o / l^hPnl"^ + o / \^hPn+l " ShPn\'^- 

Jn -^ Jn -^ Jn ^ Jn 

Moreover, the second integral on the left-hand side of (4.12) can be written in terms 
of ShPn+i in an obvious way. Finally, by splitting /' into f[ + f'2, observing that the 
contribution due to the terms involving f[ is nonnegative since /{ is monotone and moving 
the other ones to the right-hand side, we see that (4.12) yields the inequality 

- / \5hPn+l\^ - ^ / |4Pn.r + 77 / \^hPn+l - ^hPn\^ + h I \6hV pn+l\^ 

"^ Jn "^ Jn "^ Jn Jn 



<- / {f2ipn+2)- f2ipn+l))5hPn+l+ / {Pn+ig' (pn+l) - Pn.g' ipn))ShPn+l- (4.13) 

Jn Jn 

The first term on the right-hand side of (4.13) is easily treated in the following way: 

(/2(Pn,+2) - f2ipn+l))5hPn+l <ch \6hpn+l\'^- (4.14) 



On the other hand, we have 



(Pn+ig'ipn+l) - Pn9'ipn))ShPn+l 

Pn+l{g'iPn+l) - g'iPn))5hPn+l + / (Pn.+ l " Pn) g' iPn)ShPn+l 

Jn 

< Ch flri+l\ShPn\ \ShPn+l\ + / (/in+l " Pn)g' {Pn)5hPn+l ■ 



Next, we deal with the last integral by using equation (2.27). Owing to our assumptions 
on g, we obtain 

/ {pn+l - Pn)g'{Pn)ShPn+l = "^ / . , " {Pn+1 ^hln " A/i„+i) (5ftP„+i 

Jn Jn ^ + ^In 



< ch Hn+l \5hPn\ \8hPn+l\ -hi V/i„+i " V ( -^ 

Jn Jn ^^ 

< ch /i„+i \ShPn\ \ShPn+l\ 

Jn 



+ 27„ 



ShP' 



'n+l 



h I ^ , " V/i„+i ■ V6hpn+i -hi 6hpn+i V/i„+i ■ V- 



1 + 27n Jn 1 + 27n 

We treat the last three terms separately. Thanks to the Holder, Sobolev, and Young 
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inequalities, and our assumption on g, we have for every o" > 



h / /i„+l \6hPn\ |4Pn+l| < c/l ||/i„+i||4 ll^Pnlli/ \\5hPn+ 



1||4 



ch 



< ah\\Shpn+l\\v -^ ll/^n+l|lv W^hPnllH (4-15) 

a 

-h / " Vpn+l ■ ^^hPn+l <ch |V/i„+i| \V6hpn+l\ 

ch 

< ah\\WShPn+i\\H H llV/in+illn (4.16) 

a 

-h I 6hPn+l V/i„+l ■ V . ^ " < ch \ShPn+l\ |V/i„+i| |Vp„| 

Jn 1 + -^7n Jn 

< ch\\6hPn+lh W'^Pu+iWh llVpnIU 

< (rh\\6hpn+i\\v + — l|V/i„+i||H {WPuWh + IIAPnlln) • (4-17) 
Now, we rewrite (4.8) as 

II^Pn|li+ UnWn < c{\\6hPn-l\\H + HPn-lH^ + l) fom = 1, . . . , A^, 

and note that we can allow the choice n = provided that we define 

p_i := po and, e.g., p_i := 0. 
Hence, we can improve (4.17). By using (4.3) and (4.7) as well, we have 

g'ipn) 



-h / (5/jP„,+i V/i„+i ■ V^ , ^ 
Ju 1 + 27n 

ch 

< ah\\5hPn+l\\v H llVpn+llln (llPnlll- + ll^Pn-llll- + llAtn-lUn + l) 

ch 

< ah\\6hpn+i\\v + — llV/in+illn {\\^hPn-i\\H + l) • (4-18) 

By recalling all these estimates, we see that (4.13) yields 

- / |4p„+i|^-- / l^hPnl"^ + 7; / \ShPn+l - ShPnl'^ + h / \6hVpn+l\'^ 

"^ Jn '^ Jn '^ Jn Jn 

f ch 

<ch \ShPn+l\'^ + 3(^h\\6hPn+l\\v ^ l|Pn+l||y ll^^iPnll^ 

Jn o- 

ch 

+ — \\Vpn+l\\H{\\ShPn-lfH + l). 

a 

Now, just by changing the value of the constant c in front of the first integral on the right- 
hand side, we can replace the last integral on the left-hand side by ||(5/iPn+i||y. Then, we 
choose a = 1/4 and rearrange. We obtain 



|2 

\v 



<ch IShpn+ll"^ + ch\\pn+l\\v\\^hPn\\H + ch\\Vpn+l\\H {¥hPn-l\\H + l) 

Jn 
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At this point, by assuming m < N — 1, we sum over n = 0, . . . ,m — 1 and have 

m— 1 „ , m— 1 



1 /* 1 — /* I, — 

- / |4Pmr + - V / \Shpn+l - ShPn\^ + -r V] PhPn+1 ||y 

"^ -^f^ ^ n=0 -^^ ^ n=0 

^ « m~l „ m— 1 

<0 / l^/iPor + c/lV / \5hpn+l\'^ + chy2\\f^n+l\\l\\5hpn\\H 
^ J^ n=0 -^^ n=0 

m—1 m—1 

+ ch^\\Vpn+lfH¥hPn-l\\]l + ch^\\Vpn+l\\]l. (4.19) 



n=0 n=0 



The second and the last terms on the right-hand side of (4.19) have been already estimated 
by (4.6) and (4.2), respectively. To treat the first term, we write (2.28) with n = and 
add Apo to both sides. Then, we multiply the resulting equality by 5hPo and integrate 
over Vl. After a rearrangement, owing to (2.4) and the assumptions on the initial data 
(see (2.7), in particular), we obtain: 

\5hP,? + h f \V6,po\' + f (/((pi) - f[{po))SHPo 
Jn Jn 

(Apo + Pog'ipo) - /i(po) - f2ipi))^hpo < cW^hPoWn- 



As (fiipi) — fi{po))(^hPo > due to the monotonicity of /{, we immediately deduce that 

\\ShPo\\H + h\\VShpo\\jj<c. (4.20) 

In particular, the desired estimate for ||(5/,po||_ff is achieved. Therefore, on recalling that 
ShP-i = because p_i = po, we see that (4.19) yields: 



\\Shpm\\H+ ^\\5hpn+l - ShPuWh + h'^WShpn+lWv 
n=0 n=0 

m—1 m—l 

< C + ch^\\Hn+l\\v PhPnllH + ch^\\Vpn+l\\H PhPn-l\\H 

n=0 71=1 

m—1 m—2 

< C + c/l ^||p„+i||y ||5/iP„||^ + c/l ^||Vp„+2||^ PftPnll?/ 

n=0 n=0 

m—1 

<Ci+C2h J2{\\Pn+lfv + \\Pn+2\\l) WShpnfn 



n=0 



for m = 0, . . . ,N — 1. Hence, we can apply the discrete Gronwall lemma (see (2.38), 
where N is to be replaced here by A^ — 1) and deduce that 



m—1 m—1 

2 

V 

n=0 n=0 

m—1 



\\ShPm\\H + y^JShPn+1 -ShPnWn + h^\\ShPn+l 
n=0 n=0 

m—1 

< Ciexp(C2h^[\\jJ,n+l\\v + ||Pn+2||y) j 



n=0 
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for m = 0, . . . ,N — 1. Owing to (4.2), we infer that 

\\ShPm\\H + y^P/iPn+1 -ShPnllH + h^\\5hPn+i\\v < ^ for m = 0, . . . , A^ - 1; 

n.=0 n=0 

moreover, on using the estimates of 6hPo and V6hPo given by (4.20), we conclude that 

N-l N-1 

max \\Shpm\\H + y^\\ShPn-ShPn~l\\H + hy^\\Shpn\\v <C. (4.21) 

m=0,...,A''— 1 ■'— ' ■'— ' 

n=0 n=0 

In particular, (4.21) yields: 

\\dtPh\\L'^iO,T;H)nL^O,T;V) < C. (4.22) 

Fifth a priori estimate. We improve (4.10)-(4.11). Owing to (4.9), on using (4.21) 
in addition to previous estimates, we immediately obtain (cf. also (2.7)) that 

\\pm\\w + Um\\H<c ioT m = 0, . . . , N , (4.23) 

l|PfellL°°(0,T;Vl/) + IIP/iIIl°°(0,T;H') + l|P/i|lL°°(0,r;Vy) — ^' (4.24) 

as well as an estimate for, e.g., ^f^ in L°°{0,T; H). 

Sixth a priori estimate. We rewrite (2.27) in the form 

(l + 7n + ln+l)hPn " ^Pn+l = - Pn ^hln ■ 

We test this equality by {pn+i — Pn), and integrate over Q. We obtain 

^ / (l + 7n + 7n.+l) \ShPn\'^ + / (V/i„+i - V/i„) ■ V/i„,+i = -/l / /i„ Sh'Jn ShPn- 

Jn Jn Jn 

As g is nonnegative and Lipschitz continuous, we infer that 

h / \6hPn\^ + - / |V/i„+iP - - / \Vpn\^ + 2 / |V(/in+l ' Pn)\^ 
<chpn \6hPn\ \5hPn\ < c/l ||/i„||4 P/iPnIU PhPnlh 

Jn 

< 2 ll^'^'""!!^ + ch\\Shpn\\v (l|V/in||?f + WPuWh) 

< - WShPnWn + ch\\6hPn\\v W"^ Pnlln + ch\\6hPn\\v y 

the last inequality being due to (4.2). By rearranging and summing over n = 0, . . . , m — 1 
with 1 < m < A^, we get: 

, m—l _ , 9 m— 1 



5^11 Vnll?/ + 2 l|V/im,||H + y Xll|V<5;,/i„,||^ 

n=0 ra=0 

_ m—l m—l 

< 2 ll^'^oll^ + ^^X^II^^^'^ll^ l|V/in||^ + c/l^||5fepn 

n=0 n=0 

m—l 

< C + ch^\\6hPn\\v l|V/i 



m—l 
n=0 
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where we have used (4.21). Now, we first apply the discrete Gronwall lemma (2.38) and 
then account for (4.21) once more. We obtain, for ?7i = 1, . . . , iV, 

771—1 -777—1 

/l^||Vn||^+||V^„||^ + /l2 Y,\\^^hf^nfH<C. (4.25) 

77=0 77 = 

Next, by (2.27), the Holder and Sobolev inequalities and the Lipschitz continuity of g, we 
infer that 

\\AHn+l\\H < c{\\6hfin\\H + \\fJ'n+lShln\\H) < c{\\6hHn\\H + H/in+lIU P^TtiIU) 
< c[\\6hiJ^n\\H + ||yU77+l||4 II^PtiIU) < c{\\6hfJ^n\\H + ||At77+l||v ||^?7P77|| v) ; 

note that in the last product we can ignore the factor ||/i„_|_i||y, due to (4.3) and (4.25), 
provided we update the last value of c. By squaring, summing up, and multiplying by h, 
we thus obtain for ?Ti = 1, . . . , A^ the estimate 



777—1 777 — 1 771—1 

|2 
\V ' 

77 = 77=0 77=0 



/l^||A/i„+i||^ < ch^\\5hfln\\H + ch^\\ShPn\ 



and we can replace the if- norm of Afin+i by the l/T-norm of fin+i thanks to (4.2). We 
collect this and (4.25) and account for (4.21) and /io G V. We have: 



Af-l TV-l 

hy^\\6hf^n\\H+ niax ||V/i„||^ + /i V'll/in+ill^ < c, (4.26) 

•'— ' n=0,...,N •'— ' 

71=0 77=0 

SO that 

\\dtJih\\L^iO,T;H) + \\'Ph\\L^{0,T;V)nL^{0,T;W) + ||P^|U°°(0,T;y) + ||Ph|| L°°(0,T;y) < C. (4.27) 

We note that (4.25) also gives the non-sharp estimate 

h\\'^dtflh\\l2^0^T:H) < c- (4.28) 

Limit and conclusion. By standard weak compactness results, we find some conver- 
gent subsequence for the interpolants. Therefore, in principle, it is understood that the 
convergence that we refer to holds for a subsequence. However, once we prove that the 
limit we find is the solution {fi, p) to problem (2.12)-(2.14), then the whole family of in- 
terpolants is convergent, due to uniqueness. For the reader's convenience, we select some 
estimates among those we have proved in the previous steps. These are: 

\\'Ph\\L°°{0,T;V)nL'2{0,T;W) + ||/^^|| L°°(0,T;y) + \\Ph\\L°°{0,T;V) < C (4.29) 

||P/7||l°°(0,T;V1/) + ||P^||L°°{0,T;iy) + ||Ph. ||l°°(0,T;VF) < C (4.30) 

\\dtPh\\L^{0,T;H) + \\dtPh\\L°^'{0,T;H)nL'^{0,T;V) < C. (4-31) 

Now, we observe that (4.31) and (2.25) imply that 

ll/Ltfe - At/ill L2(0,T;H) + IIP^ - Ph\\L^{0,T;H) <ch (4.32) 

W'Ph - Ph\\LHO,T;V) + IIp^ - ph\\L^O,T;V) <ch. (4.33) 
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This yields, in particular, that the weak limits we find for /I^^, /x , and /i/i, by using (4.29) 
and weak compactness results coincide and that the same happens for p^j, p , and ^h- 
Therefore, we can conclude that some functions ^ and p exist such that 

71^, fi , Jlh ^ p weakly star in L°°(0,T; V), (4.34) 

/Ift -^ fi weakly in ^^(0, T; W), (4.35) 

'Ph^ Ph' Ph^ P weakly star in L°°(0, T; W), (4.36) 

dtjlh -^ dtfi weakly in L'^{0, T; H), (4.37) 

dtph -^ dtp weakly star in L°°(0, T; H) D L^{0, T; 1/) . (4.38) 

Now we prove that {p, p) satisfies (2.8)-(2.11) and solves problem (2.12)-(2.14). 

We remark that the topology alluded to in the statement of Theorem 2.4 is precisely 
the topology associated with the convergences specified in (4.34)-(4.38). Clearly, (2.8)- 

(2.10) are fulfilled. Moreover, the Cauchy conditions (2.14) are satisfied, because iph,Ph) 
converges to {p,p) at least weakly in C°([0,T]; if). Therefore, it remains to check that 

(2.11) holds and that equations (2.12)-(2.13) are satisfied. To do that, we read the discrete 
problem (2.27)-(2.28) in terms of the interpolants. We have: 

(1 + 2^JdtPH + Jihdt^h - A71, = 0, (4.39) 

dtPH-^Ph + nPk) = f£jiP,)- (4.40) 

Hence, the main problem consists in identifying correctly the limits of the nonlinear terms 
and those of the products. To this end, we recover some strong convergence (without 
looking for sharpness, since it is not necessary). We first recall that the embeddings 
V d H and W C C°(r2) are compact, so that we can apply [32, Sect. 8, Cor. 4] and 
deduce that 

Ph -^ P strongly in C°([0, T]; H) and a.e. in Q, (4.41) 

Ph^p strongly in CO([0,T];CO(n)) = CO(Q). (4.42) 

By combining this with (4.32) and (4.33), we infer that 

/i/i, p^ p and p^, p -^ p strongly in L^(0,T; H) and a.e. in Q. (4.43) 

We point out that a.e. -convergence actually holds for a subsequence. Take now p. G (0, p*) 
and p* G (p*, 1), where p^,^p* G (0, 1) are given by Theorem 2.1. Then, (4.42) implies that 

Pm ^ Ph ^ P* in Q, provided that h is small enough. (4.44) 

This means that the same bounds hold for p„, n = 0, . . . ,N (where {pn)n=o is the vector 
associated with ph), and hence also for p^ and p . As /', g, and g' are Lipschitz continuous 
on [p,,p'], we deduce a bound for ||/'(p/i)||oo and for the L°° norm of the other functions 
we are interested in. By combining this with the a.e. -convergence implied by (4.43), we 
deduce that 

0(Ph)> (f>ipf) -^ 0(P) strongly in Lp{Q) for p G [1, +00) and (j) = /', g, g'. (4.45) 

In particular, (2.11) holds. On the other hand, we also have 

\dt%\ = \5hln\ = \Shg{pn)\ < c\5hpn\ = c\dtph\ a.e. in /„+i, for n = 0, . . . , A^ - 1 , 
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so that (4.31) yields an estimate of dt^h in -Z^°°(0, T; H). Hence, thanks to (2.24), we have 

Wlh -7hlU°°(0,T;/f) < ch\\dah\\L'^{0,T;H) < ch, 

whence even ^y^ converges to g{p), e.g., strongly in L^{Q). Therefore, we infer that 

dah -^ dtgip) weakly star in L°°(0, T; H). (4.46) 

Finally, as to the limits of the products in (4.39)-(4.40), we can infer that 

Ih^tJih ^ 9{p)dtfi, Jl^dah^ fidtgip), Aih^'(P^) ^ ^^'(p), weakly in Li(Q). 
Therefore, (2.12)-(2.13) follow from (4.39)-(4.40), and the proof is complete. 



5 Proof of Theorem 2.5 

In this section, we prove Theorem 2.5. It is understood that h is as small as needed; 
oftentimes, we do not pause and quantify such smallness precisely. First of all, we remind 
the reader that the interpolants ph, p , and p^ are uniformly far for and 1 (see (4.44) 
and the subsequent lines). Therefore, without loss of generality, we can assume that the 
derivative function /' is Lipschitz continuous. We need additional a priori estimates. 

Auxiliary a priori estimates. We prepare an estimate for ||V(5/iPo||h- To this end, 
we notice that (2.28) with n = can be written as 

6hpo - hA6hpo = f'ipo) - fi.pl) - ^0 , (5.1) 

where t/'o := — ^Po + f'iPo) — Pog'iPo)- As '?/'o £ ^ by (2.33), we can test (5.1) by —A6hPo 
and integrate by parts. In view of the Lipschitz continuity of /', we find out that 

||V4po||h + h\\^ShPo\\H < ch\\6hPo\\H \\A6hPo\\H + W^Mh ||V5ftPo||i? 
< - \\A6hPo\\H + ch\\6hPo\\H + I II V^^PoIIh + c 

By accounting for (4.20), we obtain the desired estimate 

\\VShPo\\h<c. (5.2) 

Let us come now to the basic estimate we need. We write (2.28) with [n + 1) in place 
of n, and take the difference between the so-obtained equality and (2.28) itself. Then, 
we multiply this difference by —A6hPn+i and integrate over Q. We easily have, for n = 
A - 2, that 



(VShPn+l - ^ShPn) ■ V6hPn+l + / (Ap„+2 " ^Pn+l) A5hPn+l 

Jo, 

{f'iPn+2)-fiPn+l))hA5hPn+l) 
+ I ifJ'n+ig'iPn+l) - fing'{Pn)){-A6hPn+l)- (5.3) 
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By the elementary identity (2.36), the first integral is equal to 

\V6hpn+l\^ - ^ / |V(5/,p„P + - / \V6hpn+l-'^Shpn\ 

On the other hand, we obviously have that 



{Apn+2 - Apn+l) A5hpn+1 = h \A6hpn+l\ ■ 

n Jn 

Now, we deal with the right-hand side of (5.3). By Lipschitz continuity, we deduce that 

- / (/'(Pn+2) - f (Pn+l)) {-A5hPn+l) <C I \pn+2 " Pn+l\ \A5hPn+l\ 

Jn Jn 

= ch \6hpn+l\\A6hpn+l\ < 7 / \A6hpn+l\^ + C h / \Shpn+l\'^- 

Jn 4 Jq Jq 

As far as the last term of (5.3) is concerned, we combine the above elementary argument 
with the Holder and Young inequalities and the Sobolev embedding V C L^{fl). We find: 

/ {Pn+l9'{pn+l) - Pn9'{pn)){-A6hpn+l) 

Jn 

<C {pn+l\Pn+l - Pn\ + \Pn+l - Pn\)\AShPn+l\ 

Jn 

< C h[\\fin+l\\4\\ShPnh + \\ShPn\\H)\\A6hPn+l\\H 

<- / \A6hPn+l\'^ + ch\\pn+l\\v\\ShPn\\v + ch\\6hPn\\H- 

4 Jn 
By collecting the inequalities we have obtained, we see that (5.3) yields: 

- / |V(5/,p„+ip-- / |V4PnP + 7^ / |V5^p„+l-V5,,p„P + - / \A6hPn+l\^ 

'^ Jn "^ Jn '^ Jn ^ Jn 

<ch l^hPn+ll"^ + ch\\pn+l\\v\\^hPn\\v + ch\\6hlJ.n\\H- 

Jn 
At this point, we sum over n = 0, . . . , ?n — 1, with 1 < ??7, < A^ — 1, and deduce that 

m—l „ , m—1 



- / |V5/,p™|2 + - V / \VShPn+l-^ShPn\^ + - V / \A6hPn+l\ 

'^ Jn ^ „^Q Jn ^ . n J'^ 

+ c max \\pn+i\\v hy^\\6hPn\\v + chy^\\ShPn\\H- (5.4) 

n=0,...,iV— 1 •^— ' ■^— ' 



n=0 n=0 



The first term on the right-hand side of (5.4) is estimated by (5.2); all other terms on 
the right-hand side have been estimated already (cf. (4.22) and (4.27)). Therefore, by 



20 Time discretization of a nonstandard viscous Cahn-Hilliard system 
recalling also (4.21), we conclude that 

max \\Shpm.\\v + hy^\\A6hpn+i\\H <c, (5.5) 

m=0,...,A' — 1 •'— ' 

n=0 

\\dtPh\\L^{0,T;V) + \\/^dtPh\\L^{0,T;H) < C (5.6) 

Consequence. In view of the regularity theory for elliptic equations and the continuous 
embedding W C L°°{Q), we derive from (5.6) that 

ll^tP/i||L2(o,T;W) < c and \\dtPh\\L^O:T;L^in)) < c (5.7) 

Moreover, as the second (5.7) means an estimate of the difference quotients associated 
to the vector {pn)n=oy ^^'^ ^^ 9 ^^ Lipschitz continuous, a similar estimate holds for the 
vector (g{pn))n=o (^^^ (^-l)), and we infer that 

\\dt%\\L^{0,T;L^{n)) < c. (5.8) 

Furthermore, by applying (2.25), we see that (5.6) also implies that 

||A(p^-p^)lU2(o,T;H)<c/i. (5.9) 

Proof of Theorem 2.5. A possible strategy could be the following: to multiply the 
difference between (4.39) and (2.12) by (jlh — A*), and the difference between (4.40) and 
(2.13) by dt{ph — p)', then, to sum up and start estimating. However, in order to split 
calculations and give more transparence to the proof, we prefer to proceed with those 
pairs of equation separately, and collect the inequalities we obtain later on. So, we first 
consider just one couple, for instance, (4.40) and (2.13). We multiply their difference 
by dti'ph — p)i integrate over Qt, where t G (0,T) is arbitrary, and add the same integral 
to both sides for convenience. We obtain: 

fJ^\dt{PH- p)? + \\\{PH- p)mi 
-A(p,-p,)-(r(p,)-/'(p)) 



-'n 



+ 9'ip.)iPu -p) + pia'ipJ - g'ip)) + {ph - p) \dt{ph - p) 



<\ fl \dt{pH-p)\' 



'0 Jn 




'0 Jn 



l^iph - pJl' + [Ph - pf + K - P\^ + \P^ - P\^ + \ph - pI'}- (5.10) 



In the above inequality, we have used the Lipschitz continuity of /' and g', and the 
boundedness of p. Now, we estimate the last integral of (5.10). Thanks to (5.9), we 
have that 



Jn 
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On the other hand, owing to (4.33), we obtain: 



\ph - P? + \p_h - P? + \ph - pH < C / / {\ph - Ph? + |P. - Ph? + \ph - P?) 



iQJn JQJU 

t 

"2 t „ I I \^ „|2 



<ch +c / \ph- P 
Jo Jn 



Similarly, we have, by (4.32), that 

l/i^ -/i|^ < c / / {\p^-ph\'^ + \Ph -pH <ch'^ + cl I Ip/j -/i|^ 



'oJn JQJn JoJn 

By collecting the above inequalities, we see that (5.10) and the Gronwall lemma yield: 

f[\dt{PH-p)\'+\\iPk-pmrv<c{h'+ f[\PH-p\'+ f I \Ph-P?] (5.11) 

Jo Jq *- Jo Jn Jo Jn ■' 

for every t G [0,T]. Now, we deal with equations (4.39) and (2.12). For the reader's 
convenience, by recalling that 7^^ = g{ph) and 7 = g{p,) (see (4.1)), we rewrite the 
former in a different way, namely, 

(1 + 2g{p^))dti2h + pAlh - ^Ph = 0. (5.12) 

Next, we take the difference between (5.12) and (2.12) and write it as 

(1 + 2g{ph)) dtiJih - p) - ^{Ph - p) + {Ph - p) 
= -2dtfi{g{ph) - g{p)) - dah [ph - P) - pdt{^h - g{p)) 
+ '^{9{ph) - g{pf))dtPh - ^{Ph - Ph) + {Ph - p) ■ 

Finally, we multiply this equality by {jlh — p) and obtain the following identity: 

dt{{\+ g{ph)) {Ph - pf} - ^{Ph - p) {Ph - p) + {Ph - pf 
= dtg{ph) {Ph - pf - 2dtp {g{ph) - g{p)){ph - p) 

- dtlh (Ph - P) {Ph - ^)- pdt (% - g{p)) {Ph - p) 

+ '^{aiph) - 9{pJ)dtPh {Ph - p)- ^{Ph - Ph) (Ph -p) + {Ph - p)^- (5.13) 

At this point, we integrate over Qt- As g is nonnegative, we get: 

If \iPh - p){t)\' + f\\{ph-p){s)ryds<j2m^ (5.14) 

^ Jn Jo j^;^ 

with an obvious meaning of Jj(t), j = 1, . . . , 7. Now, we estimate these integrals, but the 
last one. By combining Holder, Young, and Sobolev, inequalities, and in view of (5.6), 
we have that 

hit)<c / \\dtph{s)\\4\\iPh- p){s)\\H\\iPh- p){s)\\4ds 
Jo 

<c \\dtPh{s)\\v\\{Ph-p){s)\\H\\{Ph-p){s)\\vds 
Jo 

<a WiJlh- p)is)\\vds + c^ / WiJlh- p){s)\\Hds, 
Jo Jo 
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where a > is arbitrary. Similarly, we infer that 

l2{t)<2 [ \\dtii(s)\\H\\{ph- p){s)\U\\{Jih~ f^){s)\Uds 
Jo 



Notice that, by means of the Gronwall lemma, we shall be able to control the last integral 
in terms of the L^(0,T)-norm of the function s i— )• ||c?t/i(s)|||^ (cf. (2.8)). We use a similar 
procedure for the next integral and notice that the same remark holds, due to (5.8). 
Indeed, we have that 

hit) < f \\dah{s)\U\{i2H- p){s)\\hUJI,- p){s)\\Hds 
Jo 

< / \\dtlh{s)\\oo\\{ph - p){s)\\H{\\{ph - p){s)\\h + \\{Ph - Ph){s)\\H) ds 
Jo 

< t\\daH{s)\\U\{p, - p){s)\\lds + c fiUri, - p){s)\\l + IK/I, - i2,){s)\\j,) ds 
Jo Jo 

<c [\\\daH{s)\\l + l)\\{i2,-f,){s)\\j,ds + ch\ 
Jo 

where the last inequality is due to (4.32). In order to treat I^it), we prove a preliminary 
estimate, namely, that 

\dt{lh-9ip))\ <c{\ph-ph\ + \P,^-Ph\ + \ph-p\} \dtph\ +c\dtiph-p)\ (5.15) 

a.e. in Q. As we argue pointwise, we fix {x,t) a.e. in Q and choose n such that t belongs 
to the interval {nh, {n + l)h]; in order to simplify the notation, we omit writing at what 
point (x, t) we work. By the mean value theorem, we find r between p„ and Pn+i such that 

a f^ / N\ 9{Pn+l) — g{Pn) ,, NO 1/ .Pn+l—Pn // no 

dt[lh - 9{p)) = j^ 9 [p)dtp = 9 [r) g {p)dtp 

= 9'{r)dtPh - 9\p)dtp = {gir) - g\p))dtph + g'{p) {dtPh - dtp) . 

As g' is bounded and Lipschitz continuous, we infer that 

\dt{lh - g{p))\ < c\r- p\ \dtph\ + c\dtph - dtp\. 
On the other hand, we have 

\r - p\ <\r - Pn\ + \pn - p\ < \Pn+l - Pn\ + \Pn " p\ 

= \Ph - Ph\ + \Ph ~ P\^ \Ph - Ph\ + 2|p^ - Ph\ + \ph - p\- 

Hence, (5.15) follows, and we can use it to estimate ^{t). We also account for the 
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boundedness oi fi and for identity (2.24) and the analogue identity concerning 2;^. We have: 

hit) <c {\ph-ph\ + \Ph ~ P^\ + 1^^ ~ ^1} \^^P^\ \^h ^ ^1 

Jo Jn 

+ c / / \dt{ph - p)\ \^h -/i| 



Jo Jn 

< c fiWMs) - MsWh + IIP,(^) - MsWh + WMs) - Pis)\\l} \\dtMs)\\l ds 

Jo 

+ 2/ \\9t{ph - p)is)fH ds + c \\(jlh-fi)is)\\jjds 

< ch'^\\dtPh\\l^{o,T;H) \\9tPh\\l^o,T;L--{n)) + c / ||5tPh(s)||L \\ph{s) ' p{s)\\h ds 

Jo 



t rt 



+ II \\dt{ph - p){sWh ds + c I \\(jlh-fi)is)\\jjds 



furthermore, estimates (4.31) and (5.7) allow us to infer that 

l4it)<ch^ + cJ\\dtMs)\L\\Ms)-pis)\\lds + ^J\\dtiph-p){s)\\lds 

+c[\i + \\dM^)\L)\m.-^^){s)\\ids. 

Jo 
Next, we deal with I^it). By accounting for (2.25) and (5.7), we deduce that 

h{t) <c \\{ph - pJis)\\^\\dtflh{s)\\H\\{Jih - f^){s)\\H ds 
Jo 

< ch^\\dtPh\\h{o,T;L--{n)) + c / \\dtJlh{s)\\jj\\{Jlh-fi){s)fHds 

Jo 



<ch +c \\dtPh{s)\\H\\{fih- l^){s)\\Hds. 
Jo 

We note at once that we shall be able to control even the last terms of the last two 
estimates with the help of the Gronwall lemma, in view of (5.7) and (4.31), respectively. 
Finally, thanks to (2.25) once more, we have: 



h{t)= [ /v(/2,.-7lJ-V(/2;, 
Jo Jn 



fi) 



Jo Jn Jo Jn 

< CT / / |V(/2/, -/i)p + c^/i^||V(9i/I/,||i2(o,r;H) 
Jo Jn 

<(T / /|V(/2/.-/i)|' + c,/i, (5.16) 

Jo Jn 

where the last inequality is a consequence of the non-sharp estimate (4.28). We stress 
that Iq is the only term of order h instead of h"^. At this point, we collect all the estimates 
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of the integrals Ij we have obtained, and come back to (5.13)-(5.14). If we choose a small 
enough, we conclude that 

1 1 /■* 



t 



<c\h+ / {i + \\dMs)\\lo + \\dtMs)rH)\\{rik-^^){s)\\lds 

\dt{pH - p)\'} + I j\dt{pH - P){s)rn ds (5.17) 





i 



for every t G [0,T]. Now, we revert to (5.11) and add it to (5.17). After rearranging, we 
apply the Gronwall lemma and obtain (2.30). This concludes the proof. 



6 Proof of Theorem 2.6 

As is clear from the proof of Theorem 2.5, to obtain estimate (2.32) there is just one 
step to modify, namely, the estimate of Jg (see (5.16)), which was based on the non-sharp 
inequality (4.28). Thus, we only have to prove that our further assumption (2.31) imphes 
that Iq must be of order h'^, not h. Moreover, it is clear that this is true whenever we 
improve (4.28) and replace it by 

N-l 

||Vc}i/i/j||^2(o,r;ii-) < c, i.e., h^\\VShPn\\H < c. (6.1) 

n=0 

Hence, it suffices to prove (6.1). In order to make our argument transparent, we prove 
some additional estimates, the first of which holds under assumption (2.31). 

Further a priori estimates. We prepare an estimate of ||5/iyUo||H. In view of (2.31), 
we write equation (2.27), with n = 0, in the form: 

(1 + 2-fo)6hPo - hA6hPo = A/io - Pi Shlo, 
and test it by 6hPo- As 70 is nonnegative, we immediately arrive at 

\ShPo\'^ + h \V6hPo\'^ < {\\Aiio\\h + ll/xilU llf^ftPolk) H^/x/^olk • 
ci Ja 

Thanks to (2.31), the Sobolev inequality, (4.26), and (5.5), we deduce that 

\\ShPo\\H + h [ \V5hPo\^<c. (6.2) 

Jn 

Let us come to the basic estimate we need. We improve (4.21) and obtain a bound for 
the second- difference quotients 6lpn (see (2.16)). We write (2.28), with (n+ 1) in place of 



Colli — Gilardi — Krejci — Podio-Guidugli — Sprekels 25 

n, and test the difference between the resulting relation and (2.28) itself by {ShPn+i—^hPn)- 
We find: 

/ \ShPn+l - Shpnl"^ + h V5hPn+l ' ^{^hPn+l " hPn) 

Jn Jn 

= - (/'(Pn+2) - f'{pn+l)){ShPn+l - Shpn) 

Jn 

+ h {g'iPn+l)5hPn + Pn5h{g' {Pn))){5hPn+l - ^hPn) 

Jn 

<Ch {\ShPn+l\ + \5hPn\ + \Pn\ \ShPn\)\5hPn+l - 5hPn\ ■ (6.3) 

Jn 
By the elementary identity (2.36), we have: 

VShPn+l ■ ^{hPn+l - ShPn) = -\V6hPn+l\^ " ^I'^hpnl'^ + ^{"^ShPn+l " V(5/iP„p. 

On the other hand, by the Sobolev inequality, (4.26), and (5.5), we infer that 

Ch {\6hPn+l\ + \5hPn\ + \Pn\ \5hPn\)\ShPn+l - ShPn\ 



'n 

< O / \^hPn+l - ShPnl^ + ch^{\\6hpn+l\\H + W^hPnWH + \\Pn\\l\\ShPn\\l) 

^ Jn 

< - / \^hPn+\ -<^hPnr + c/l^(||(5ftP„+i||^+ ^^hPufk + ||/in.||y H^hPnUv^) 

2 Jn 

^ Jn 

Now, we combine this estimate, the identity just above, and (6.3). Then, we divide by h 
and sum over n = 0, . . . , m — 1, where 1 < m < iV — 1. We obtain: 



h 



m—\ ^ -, 771—1 



I I h J_ J -.(ft J. 

^mPnfn + 2 l|V4Pm||^ + 2 ^\\^^hpn+\ " ^ ^hpufu 

n=0 n=0 

_. m—l m— 1 



n=0 n=0 

At this point, by (5.5), (4.6), and (4.26), we conclude that 

N-2 

hJ^W^lpufH < c. (6.4) 

n=0 

Consequence. With a view toward deriving an estimate for (5^7n, we begin by arguing 
pointwise. So, for a. a. (x, t) E Q (once again we omit writing at what point of Q we work) 
and for suitable ri between pn+2 and p„+i, and r2 between p„ and Pn+i, we have by the 
Taylor formula: 

I^h7n| = h-^ \9{pn+2) - g{pn+l) + g{pn) " 9{pn+l)\ 

= h~^ \g\pn+l){.Pn+2 - Pn+l) + \ fi'"(^l)(Pn+2 " Pn+lf 

+ fi''(Pn+l)(Pn, - Pn+l) + |fl'"('"2)(Pn " Pn+lf\ 
< C\5lpn\ +c{\5hpn+l? + \5hPn?)- 



26 Time discretization of a nonstandard viscous Cahn-Hilliard system 
Now, we square this pointwise estimate, integrate over Q, sum over n, and deduce that 

7V-2 N-2 N-1 

hY.\\shnrH<chj2¥iPnrH+chY,\\s,pX- 

n=0 n=0 n=0 

Then, (6.4), the Sobolev inequahty, and (5.5) yield: 

N-2 

hJ2\\^hnrH<c. (6.5) 

n=0 

Proof of Theorem 2.6. As said before, it suffices to prove (6.1). We reason that, in 
order to obtain the analogous estimate for the solution to the continuous problem, one 
first differentiates (2.12) with respect to time and then tests the resulting equality by dtfi] 

this yields the desired term / / |V9i/ip on the left-hand side. The idea is to perform 

Jo Jn 
the corresponding procedure on the discrete equation (2.27). However, it turns out that 

the calculation in the discrete case becomes simpler if one tests by the analogue of the 

product (1 + 2g{p))dtfi. To simplify the notation, we introduce the vector it defined by 

TTn ■■= (1 + 2'-fn)ShfJ'n ioi U = 0, . . . , N - 1. (6.6) 

We write (2.27) with {n + 1) in place of n, and take the difference between the resulting 
equality and (2.27) itself. Then, we test this difference by iin+i and integrate over Q. By 
taking the elementary identity (2.36) into account, we obtain for n = 0, . . . , A^ — 2 that 

T: / kn+ll^-TT / knP + 77 / kn+l - 7r„P + /?, / V5/iyU„+i ■ V7r„+i 

^ 7a "^ Jn "^ Jn Jn 



{fJ'n+2 Shln+1 - fJ^n+1 ^hln) T^n+1 ■ 

By computing the fourth term on the left-hand side with the help of (6.6), recalling that 
g is nonnegative, and rearranging, we deduce that 

/ |7r„+ip- / |7r„P+ / llTn+l - TTnl'^ + 2h / \V6hfin+l\'^ 

Jn Jn Jn Jn 

< -2 / (/in+2 Sh'-Tn+l - yUn+l 5/i7n) T^n+l - 4:h 6hfin+l V6hfin+l " V7„+l 

Jn Jn 

= -2 / nn+l{Hn+2 - fJ'n+l) 5hln+l " 2 / 7r„+i /X„+i ((5h7n+l " <5h7n) 

Jn Jn 

-Ah 5hPn+i VShPn+i ■ V7„+i 



2h / TTn+l ShPn+l Sh1n+1 - 2h TTn+l /in+1 5^7' 



-Ah 5hPn+i "^ShPn+i ■ V7„+i . (6.7) 

Jn 

Now, we estimate each term of the right-hand side separately, before summing over n, in 
order to simplify the notation. For the ffist one, we use Holder and Sobolev inequalities. 



Colli — Gilardi — Krejci — Podio-Guidugli — Sprekels 27 



and estimates (5.5) and (6.4). We have: 



2h / TTn+lShl^n+lShln+l 



< 2/l||7r,„+i||4 ||(5ft/i„+i||// ||<5h7„+i||4 < Ch\\5hl^n+l\\v Ph/^n+llln PhPn+lIk 
- 7 {W'^^hPn+lWH + ¥hf^n+l\\H) + C /l||5/,/i„+i ||| 

< - \\V5hftn+l\\H + c/i Ikn+lll^ 

< - \\V6hfin+i\\H + ch\\nn\\H + C'/l||7r„+i -TTuWh- 
For h small enough, namely, for h < 1/(3C), we conclude that 

/• 7 -1 

-2h / 7r„+i(5/i/i„+i(5/i7„+i < - ||V5h/i„+i||^ + c/i||7r„||^ + - ||7r„+i - 7r„||^. 

Next, by (4.26), we similarly have: 

~2h / 7r„+i/i„+i(5^7„ < c/i||(5;j/i„+i||4 ||yu„+i||4 ||5^7„||h 
Jn 

< -{\\V6hHn+l\\H + \\6hfin+l\\H) + C /l ||/i„+i ||y ||5^7„||^ 

< - \\V6hfin+i\\H + h i|vr„,||^ + - ||vr„+i - 7r„,||^ + ch ||5^7n||^ , 

for sufficiently small h. Finally, by accounting for (4.23), Sobolev inequality, and the 
compactness inequality (2.35), we have that, for h small enough, 

-Ah / 5/,yU„+l V5/,yU„+l ■ V-fn+l < c/i||5/,yU„+l||4||V(5/,/i„+l||H ||Vp„+i||4 

Jn 

< -\\VShlJ^n+l\\H + ch\\Shf^n+l\\l\\'^ Pn+l\\v < T II V5/,/i„+i ||^ + C /l||(5/,/i„+i II4 

< ^ \\V6hl^n+l\\H + ^(i ||V4/in+l||H + c\\6hfin+l\\H) 

At this point, we combine the inequalities just obtained with (6.7) and note that the 
terms involving vr„+i — vr^ cancel out. Then, we sum over n = 0, . . . , (rra — 1), with 
1 < m < (A^ - 1). We obtain: 

m—l m—1 m—1 

hmWn + h^\\V6h^^n+l\\H < IMh + ^ X^lknll^ + h^\\6l-fn\\H 
n=0 n=0 n=0 

and the discrete Gronwall lemma allows us to deduce that 

m-l N~2 

W^^Wh + h^\\VShf^n+i\\H < c(^||7ro||^ + c/i^||5^7„||^ 

n=0 n=0 
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for 1 < 772 < (A^ — 1). From (6.2) and (6.5), we infer that 

Af-2 



n=0 



This and (6.2) yield (6.1), and the proof is complete. 

Remark 6.1. As a consequence of estimates (6.1) and (6.4), the solution to the continuous 
problem enjoys the following additional regularity properties: 

VdtjjeL^iQ) and dfpeL^{Q). 

This can give even more: for instance, equation (2.13) can be differentiated with respect 
to time, to show that AS^p belongs to L'^{Q) as well, so as to conclude that 

peH\0,T;H)r]H\0,T;W). 

However, this regularity result could be proved formally and directly for the continuous 
problem. 
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